Pseudorapidity shape of elliptic flow as signature for fast equilibration in relativistic 
heavy-ion collisions at energies up to ^/snn = 200 GeV 
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The implications of parton recombination processes on the dynamics of ultrarelativistic heavy-ion 
reactions are investigated. To do so, the quark-gluon string transport model has been extended 
for partonic recombination and fusion processes. Parton recombination leads to short equilibration 
times and improves significantly on the theoretical description of measured directed and elliptic 
flow, i.e., Vi and V2, distributions in Au+Au collisions at ^/snn = 200 GeV, in particular what 
concerns their pseudorapidity dependence. The shape of V2{if) is found to be closely related to fast 
thermalization. 
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I. INTRODUCTION 

Among the various experimental studies of ultrarel- 
ativistic heavy-ion collisions operates the Relativistic 
Heavy Ion Collider (RHIC), at Brookhaven National 
Laboratory since 2000, to investigate gold-on-gold col- 
lisions up to y/sNN = 200 GeV. After many years of op- 
eration strong experimental evidence has been accumu- 
lated that at RHIC energies indeed a new state of matter 
is created that is qualitatively different from a hadron gas 
(see Ref. [l| and references therein). This state seems, 
however, not to behave like a weakly interacting gas of 
deconfmed partons - as could have been naively expected 
- but rather like a strongly coupled quark-gluon plasma 
(sQGP). One argument toward such a scenario is the 
large elliptic flow observed at RHIC @, |, H H H- The 
development of strong elliptic flow requires short equi- 
libration times and large pressure gradients to drive the 
dynamics, both being characteristic features of a strongly 
interacting system 0, HI • 

Another evidence for this picture is provided by the 
hadron species dependence of the elliptic flow [1, Q . The 
observed scaling with the number of constituent quarks 
can naturally be explained by the assumption that the 
elliptic flow is to most extent already created in the par- 
tonic phase and transferred to the hadrons th roug h a 
partonic recombination mechanism [Tol . [Til [T2I ll& ITil. 

Em, da. 

In the present work we demonstrate that not only the 
size but also the shape of the observed pseudorapidity 
profiles of the elliptic flow vi provide strong evidence for 
a partonic recombination mechanism in combination with 
fast thermal equilibration inside a sQGP phase. As basis 
for these investigations serves a microscopic quark-gluon 
string model (QGSM) that has been extended to allow for 
a partonic recombination procedure motivated by parton 
coalescence models [HI, 03, EH, E3, i.e., to model effec- 
tively the dynamics of a strongly coupled quark plasma 
formed in the very dense stages of ultrarelativistic heavy- 
ion collision. By construction, string cascade models do 



not contain an explicit quark-hadron phase transition. 
However, during the temporal evolution of a heavy-ion re- 
action a dense and strongly interacting plasma is formed 
within such approaches as well. The system consists of 
partons and color-flux tubes (or strings). Thus, such 
models can serve as starting base to study the dynamics 
of the dense and strongly interacting medium created at 
RHIC, in particular because transverse as well as ellip- 
tic flow at SPS energies are well reproduced within the 
string- cascade approach [l8l [lli I20L l2l1| . 

It was shown in Refs. [22I , |23[ that the standard ver- 
sion of the microscopic quark-gluon string cascade model 
(QGSM) is able to describe fairly well the bulk proper- 
ties of the elliptic flow V2 measured in ^snn = 200 GeV 
Au+Au reactions at RHIC. However, a difference in 
shape in particular between the pseudorapidity distribu- 
tions of the elliptic flow at midrapidity determined in the 
experiment and in the simulation has rankled. The shape 
of the V2 distributions have been found to be closely re- 
lated to anisotropies in the corresponding energy density 
profiles and the degree of equilibration [24[. Although 
the system reaches quickly some sort of a pre-equilibrium 
stage within standard QGSM the processes included are 
not sufficient to drive the system to a completely ther- 
malized state. As shown in the following, this changes 
dramatically when partonic recombination processes are 
included. The final consequences for the anisotropic flow 
turned out to be even more remarkable. 



II. QUARK-GLUON STRING MODEL WITH 
PARTON RECOMBINATION 

The standard Monte Carlo version of the quark- 
gluon string model (2|| that serves as the basis for the 
present investigations has been described in detail in 
Refs. [22|, [2||. Based on Gribov-Regge Theory, the pro- 
duction of hadrons is described by excitation and decay 
of open strings with different quarks or diquarks on their 
ends. In that sense the model already incorporates the 



2 



partonic structure of hadrons and therefore can provide 
a framework for the inclusion of partonic recombination 
processes. To trigger these processes, a certain critical 
(energy) density is needed to allow the hadrons to over- 
lap. Concerning the implementation of partonic recom- 
bination processes we apply essentially the method pro- 
posed in Ref. 13]: The partons of three hadrons are 
allowed to enter into a recombination process if their 
spatial distance, in the center-of-mass (cm.) frame of 
the corresponding hadrons, is less than A^ = 0.85 fm. 
From all possible triplets of hadrons satisfying these spa- 
tial constraints one is randomly chosen according to usual 
Monte Carlo methods and the corresponding hadrons are 
then decomposed into their constituent partons (quark- 
antiquark or quark-diquark). Each parton is given a mo- 
mentum fraction z of the initial momentum of its hadron 
and, additionally, the partons of each pair obtain a trans- 
verse momentum pt of opposite sign. Both, z and px are 
generated from the standard parton distribution func- 
tions used within QGSM. The (random) selection of a 
triplet of hadrons by the coordinate space criterion 
is based on the approximation that distances of hadrons 
are equal to (averaged) distances of their constituent par- 
tons. In that sense the selection procedure really acts on 
the partonic level. 

The recombination process itself requires also an over- 
lap of the participating partons in momentum space. The 
distance in momentum space between the pairs of partons 
is evaluated in the center-of-mass frame. The probabil- 
ity for recombination of two partons is then given by the 
covariant distribution II 311 
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Additionally, if the partons of three mesons participate, 
recombination of three quarks and three antiquarks is 
possible with a probability distribution fl3| 
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leading effectively to a fusion of three mesons into a 
baryon-antibaryon pair. We want to stress at this point 
that our fusion process, which allows pseudoscalar and 
vector mesons to "rearrange" to BB on the partonic level 
(no cross sections have been introduced), is qualitatively 
consistent with the method applied in Ref. [271 ] . where a 
quark rearrangement model for BB annihilation to three 
mesons (pseudoscalar and vector, which later decay to 
additional pions) has been employed. For example: The 
fusion of two p mesons and a pion into pp is kinemati- 
cally favored compared to the fusion of three pions (phase 



space). However, a more quantitative investigation has 
to be delayed to future work. 

The momentum coalescence radius A p , which enters 
in the distributions |T]) and |2j), is given by the uncer- 
tainty principle A^Ap = 1, where the actual distance 
between the participating partons is used for A^. To 
exclude recombination of particles with highly different 
momenta high-momentum cutoffs of A p max = mi + ?ti 2 
for recombination processes 3H — > 3H {H = M, B, B) 
and A Pjmax — nix + m% + + 0.2 GeV for the afore- 
mentioned fusion process 3M — * BB are used, where 
rrii corresponds to typical constituent masses. We want 
to stress that recombination and fusion happen on the 
partonic level, i.e., the notation denotes the contained 
partons, which participate in those processes, in the fol- 
lowing way: M = {qq) M , B = {qqq) B , and B = {qqq)s, 
respectively. 

The recombination (or fusion) processes can take place 
if the quantum numbers allow the partons of at least two 
of the three overlapping hadrons to recombine into new 
hadronic correlations, whereas this process is mediated 
by the partons of the third hadron. Effectively, this is a 
three-body interaction or at most a kind of in-medium 
two-body interaction but no additional two-body inter- 
action in the vacuum. Hence the total vacuum hadron- 
hadron cross section implemented in the QGSM is not 
changed. Anyway, as mentioned above no explicit in- 
medium cross sections have been introduced because all 
recombination processes are dynamically treated on the 
partonic level. Nevertheless it is possible to estimate ef- 
fective cross sections for the implemented partonic re- 
combinations between the hadronic correlations. The 
maximum value of the inclusive cross section for recom- 
bination processes like 3M — > 3M or 3M — » BB is ba- 
sically determined by the aformentioned spatial coales- 
cence radius A^. The allowed maximum value of 0.85 fm 
[HI is geometrically related to a cross section not larger 
than about 23 mb. Because the dense medium produced 
in the overlap zone of a highly energetic Au+Au collision 
is dominated by "pionic correlations" , the recombina- 
tion of partons coming from and going to two "pions" in 
presence of the partons of a third (medium) "pion" hap- 
pens most frequently. Nearly 90% of all recombination 
processes with three pionic states in the initial channel 
are reactions of this kind, whereas the "fusion process" 
3ir — > pp is distinctly suppressed. Accordingly, the effec- 
tive cross section of the latter process is several orders of 
magnitude smaller than the estimated cross section of the 
first process, which turns out to be not larger than 20 mb. 
Furthermore, we have estimated effective cross sections 
for the partons (quarks and antiquarks) involved in re- 
combinations by their distance in momentum space us- 
ing a q q ex {A p (q, q))~ 2 ■ Depending on the center-of-mass 
energy of the (anti)quarks, the partonic cross section is 
in the order of a few millibarn but definitely not larger 
than the aforementioned 23 mb. The rate of recombina- 
tion processes is determined by the cutoff parameter for 
; : increasing A p ^ max increases the 



high momenta A p 
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recombination rate. However, the additional recombina- 
tion processes with highly different momenta have only 
very small cross sections. Because the bulk of recombina- 
tion processes happens to take place with A p (q,q) close 
to the cutoff A p ^ max , an averaged cross section would 
yield a value comparable to parton scattering cross sec- 
tions used, e.g., in the AMPT model [2^ |. 

The implemented procedure always ensures that only 
physical particles can be final hadronic states, e.g., no 
final diquark states are possible. Also the reaction 
HBB — > H3M, i.e., the partons of a hadron and a 
baryon-antibaryon pair form finally a hadron and three 
mesons, which contains as a subprocess the backreaction 
to "meson fusion" , is not included in the recombination 
scheme because annihilation of baryons and antibaryons 
is already implemented in the QGSM as a standard two- 
particle reaction B + B — ► X . Furthermore, the model 
allows for the formation of resonance states and their 
decay. Consequently, subprocesses such as 2M — > 3M, 
which can also increase the particle number, are incor- 
porated in the model. However, if the outcome of a re- 
combination process would yield any unphysical particle 
it cannot take place. Either an other valid recombination 
process is then chosen or, if none is possible, no recom- 
bination occurs. 

Apart from recombining to new mesons, (anti)baryons 
or baryon-antibaryon pairs, quark- antiquark annihilation 
is possible when the partons of three mesons participate 
in a recombination process. Thus a quark-antiquark pair 
of the same flavour, but belonging to different mesons, 
may annihilate with a redistribution of its energy and 
momentum to the new mesons formed by the remain- 
ing partons. The probability of this annihilation process 
3M -> 2M with respect to 3M -> 3M, i.e., P a = 0.04, 
has been adjusted to reproduce the experimental dN/dr) 
charged hadron multiplicities. By means of this anni- 
hilation process 3M — ► 2M an effective backward reac- 
tion for diffractive scattering is included. From the pos- 
sible recombination and annihilation processes, includ- 
ing the case that nothing happens at all, i.e., all par- 
tons recombine to the original hadrons, the actual reac- 
tion is randomly chosen. These processes are checked for 
all combinations of overlapping hadrons, making thereby 
sure, however, that the particular partons of the selected 
hadron triplet can only once per time step - which is 
about 10 -5 . . . 10 -3 fm/c - participate in such processes. 
If no recombination (annihilation) or fusion processes 
take place normal elastic or inelastic scattering occurs. 

Charge conservation is automatically guaranteed 
within our approach, but conserving energy and three- 
momentum simultaneously is not possible for the re- 
combination processes described above. Generally, there 
are problems to conserve simultaneously energy and mo- 
mentum within partonic coalescence/recombination ap- 



proaches, see, e.g., Ref. [17| . where conservation of mo- 



mentum has been chosen whereas conservation of en- 
ergy has been violated. In our model, conservation of 
energy is violated in most recombination processes on 



the level of few percentages only. However, applying 
a rescaling procedure for the momenta in the center- 
of-mass system of the produced hadronic correlations, 
we are able to ensure simultaneous conservation of en- 
ergy and three-momentum with a precision better than 1 
MeV. To conserve the initial energy the three-momentum 
components of all particles involved in a particular reac- 
tion are scaled by a constant factor which is determined 
iteratively. Moreover, by such a procedure ratios of mo- 
menta, e.g., anisotropic flow coefficients, are remain un- 
changed. In other words, the anisotropic flow is not ar- 
tificially influenced by the momentum rescaling. 

All particles are allowed to interact via the recombina- 
tion procedure described above, even non-formed (pre-) 
hadrons. In this spirit the model does not create a "sys- 
tem of free partons" but effectively emulates a medium 
of very strongly correlated partons, i.e., quark-antiquark 
and (anti)quark-(anti)diquark states. To not shorten or 
modify the formation times of the nonformed "hadrons" , 
their formation time is now interpreted on the quark 
level, meaning the time for the constituent partons to 
be fully created. As a consequence of that each hadron 
gets a formation time for each constituent [quark, an- 
tiquark or (anti-) diquark] t\ and t f 2 . Within a normal 
production process such as string-breaking t{ and t\ are 
equal. For the recombination mechanism the newly pro- 
duced hadronic states get the formation times of their 
constituents. The new hadron will be fully formed if the 
lifetime exceeds the larger formation time and partially 
formed (like a leading hadron) if the lifetime exceeds the 
smaller formation time. These partially formed hadrons 
are allowed to rescatter assuming additive constituent 
cross sections. 



III. EQUILIBRATION AND ANISOTROPIC 
FLOW OF CHARGED HADRONS 



First we examine the influence of parton recombina- 
tion on the kinetic equilibration in the central cell of the 
overlap zone of Au+Au collisions. This cell is given by 
—2 fm < x,y,z < 2 fm for a central Au+Au reaction 
with impact parameter b = fm and 2 fm < x < 6 fm, 
—2 fm < y, z < 2 fm for a semiperipheral reaction 
(6 = 8 fm) , respectively. The corresponding equilibration 
ratio i?LE = (P-x + Py)/2P z is determined by the pressure 
components in x, y and z directions, P x ,y,z (see discussion 
in Ref. [IJ]). The QGSM results with and without imple- 
mentation of parton recombination are depicted in Fig. 
[T] A difference in the degree of local kinetic equilibration 
is clearly seen when one compares the standard QGSM 
results with the parton recombination scenario. The lat- 
ter mechanism leads very effectively to a much faster and 
smoother equilibration for both, the central and semipe- 
ripheral collisions. Essentially, the decomposition and re- 
combination procedure reshuffles rapidly the momenta of 
partons and finally of the hadrons which drive the system 
undoubtedly into kinetic equilibrium within short times. 
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FIG. 1: Influence of parton recombination on the kinetic 
equilibration time. The figure shows the time dependence 
of the equilibration ratio Rle in the central cell of the over- 
lap zone for central and semiperipheral Au+Au reactions at 
^/sjvjv = 200 GeV for QGSM simulations with (filled sym- 
bols) and without (open symbols) partonic recombination 
processes. 



QGSM Figure [2] shows the pseudorapidity depen- 

dence of the elliptic flow v-i of charged hadrons for min- 
imum bias Au+Au collisions at ^/s^n = 200 GeV. The 
QGSM results including parton recombination and the 
standard approach without recombination [22J are com- 
pared with the experimental data of the PHOBOS Col- 
laboration [f| . The elliptic flow obtained within the stan- 
dard QGSM displays a strong in-plane alignment in ac- 
cordance with the experimental findings. At midrapdity 
\t]\ < 1 the flow parameter V2 is almost constant, but then 
it rises up slightly followed by a rapid drop at \rf\ > 2. 
The emergence of this peculiar double bump structure in 
v-zirj) was strongly connected with the model dynamics 
(for more details, see Ref. 22]). In contrast, the ex- 
perimentally observed elliptic flow shows a pronounced 
peak at midrapidity and a steady decrease for \rf\ > 1 
0, [3Xj] . This behaviour is remarkably well reproduced 
when partonic recombination processes are taken into ac- 
count. The parton rearrangement processes in the dense 
medium lead to a redistribution of the elliptic flow of 
the final hadrons toward midrapidity, i.e., vi is accumu- 
lated at \r}\ « 0, whereas it is distinctly reduced in the 
region around \r}\ « 2. Thus, the double bump structure 
in V2(ri) obtained in the standard version of the QGSM 
disappears. 

This striking feature holds also for the centrality de- 
pendence of v-2{?l) as depicted in Fig. [3] Here, the results 
from the QGSM simulations and the PHOBOS analysis 



For central reactions, the matter in the central cell is 
practically fully equilibrated, i.e., -Rle ~ 1, after very 
few fm/c time of evolution. Here, of course, the effect 
of the recombination mechanism is strongest, because it 
is clearly a density-dependent mechanism. But even for 
semiperipheral collisions the local pre-equilibrium stage 
identified by i?LE > 1 is remarkably shortened compared 
to the scenario without recombination of partons in the 
dense stages of the expanding medium. It should be no- 
ticed that the offset of t cvo \ seen in Fig. Q] is due to the 
chosen cell size and is of technical nature. In the first 
few time steps the pressure is dominated by the longitu- 
dinal flow of the nuclei penetrating the cell. Disregarding 
this offset one can see that parton recombination reduces 
the equilibration time of the system approximately by a 
factor of 5 from ~ 10 fm/c to ~ 2 fm/c. 

The interplay between fast thermal equilibration and 
the amount of anisotropic flow, in particular the ellip- 
tic flow, measured in high-energy heavy-ion reactions is 
strongly debated (see e.g. [EJE 0, H, [1, H, HiJ [3(| and 
references therein). In Ref. [24[ the idea has been sup- 
ported that a fast and complete thermal equilibration is 
not strictly necessary to produce iarge elliptic flow. Nev- 
ertheless, a conspicuous difference in shape in particular 
between the pseudorapidity distributions of the elliptic 
flow at midrapidity determined in the experiment and in 
the simulation was observed. Hence it is very natural to 
study the effect of the quark recombination mechanism 
on the azimuthal anisotropy parameter v-i within the 
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FIG. 2: Influence of parton recombination on the pseudora- 
pidity dependence of the elliptic flow. Results for V2(rj) of 
charged hadrons from the standard QGSM (open symbols) 
[22| and those obtained with parton recombination (filled 
symbols) are compared to PHOBOS data of minimum bias 
Au+Au collisions at ^/snn = 200 GeV The statistical 
error bars are shown together with the systematic errors of 
the experimental data (gray boxes). 
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(combined data from the hit- and track-based methods) 
for three different centrality classes are overlaid. The 
QGSM, including partonic recombination processes, is 
able to describe the magnitude as well as the pseudo- 
rapidity dependence of Vz(rf) remarkably well for all the 
centrality classes, ranging from central via midcentral to 
peripheral in accordance with the definitions in Ref . [B[ . 
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FIG. 3: Pseudorapidity distributions of V2 for charged 
hadrons from ^/sjvjv = 200 GeV Au+Au reactions for three 
centrality classes according to the PHOBOS analysis The 
identification marking is the same as in Fig. [2] 

This is a highly non-trivial result, since a simultaneous 
description of both observables has neither been achieved 
by other standard string-cascade transport models such 
as relativistic quantum molecular dynamics (RQMD) or 
ultrarelativistic quantum molecular dynamics (UrQMD) 
32 , 13311 , nor by purely hydrodynamical calculations 
3J, [35|. So far, only a hydrodynamics+cascade hybrid 
approach with Glauber model initial conditions was able 
to give a fair description of the experimental data, with 
the exception of the midrapidity region in the most cen- 
tral collision class [36]. There it has been argued that 
the hadronic cascade provides the right amount of dis- 
sipation to bring the ideal fluid prediction down to the 
measured values, especially in peripheral collisions and 
away from midrapidity. 

Following this argument, the vi results obtained with 
the extended QGSM may be interpreted in a comple- 
mentary way as follows: In contrast to a highly dissi- 
pative hadronic medium, the parton recombination pro- 
cesses lead to a reduction of the mean free path in the 
very dense stages of a heavy-ion collision at midrapidity. 
Accordingly, the viscosity of this strongly interacting par- 
tonic medium is effectively lowered in comparison to the 



pure hadronic medium, i.e., the rearrangement processes 
on the partonic level reduce the amount of dissipation 
in the highly dense matter and enhance the elliptic flow, 
especially in the midrapidity region, to bring the theo- 
retical predictions in line with the data. Thus far the 
QGSM upgraded by the locally density-dependent par- 
ton recombination mechanism quasi models the possible 
dynamics of a sQGP from a microscopical point of view. 

A complementary observable regarding anisotropic 
flow phenomena is the directed flow v\ of produced par- 
ticles. In Fig. QJ the results for the pseudorapidity de- 
pendence of this flow component, i.e., v\(rf) of charged 
final hadrons, obtained by QGSM simulations with and 
without parton recombination processes are compared to 
the corresponding experimental data from the PHOBOS 
Collaboration for 0% to 40% central Au+Au collisions 
at the highest RHIC energy of ^JJWn = 200 GeV. The 
result of the standard QGSM shows a characteriatic wig- 
gle structure extensively discussed in Ref. [22|, whereas 
the scenario with parton recombination is able to repro- 
duce the experimentally observed directed flow very well 
in the broad midrapidity region where v\ is essentially 
flat and close to zero. From the microscopical point of 
view, this is another hint for a strongly interacting par- 
tonic medium which exists during the early dense stage 
of ultrarelativistic heavy-ion collisions around midrapid- 
ity. Thus, our findings are in line with the results of 
an anisotropic flow study using a multiphase transport 
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FIG. 4: Influence of parton recombination processes on the 
pseudorapidity shape of the directed flow Vi for charged 
hadrons. The results obtained with the standard QGSM 
(open symbols) [2^ ] and those from a simulation using the 
QGSM extended by parton recombination (filled symbols) are 
shown in comparison to PHOBOS data [33]. The systematic 
errors of the experimental data are indicated by gray boxes 
and the statistical errors are indicated by bars. 
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(AMPT) model that includes both initial partonic and 
final hadronic interactions [28[. There also the conclu- 
sion has been drawn that the matter produced during 
the early stage of Au+Au collisions at ^snn = 200 GeV 
in the pseudorapidity region |ry| < 3 is dominated by par- 
tons. 



IV. SUMMARY AND CONCLUSIONS 

In summary, it has been demonstrated that partonic 
fusion and recombination processes which occur in the 
very dense medium created in ultrarelativistic heavy-ion 
reactions during the early stages lead to short relax- 
ation times and drive the system to fast kinetic equi- 
librium. The basis for these investigations was a micro- 
scopic transport model, namely the quark-gluon string 
model (QGSM) based on the color exchange mechanism 
for string formation, which has been extended by a lo- 
cally density-dependent partonic recombination proce- 
dure to model effectively the dynamics of a strongly cou- 
pled quark plasma and final hadronic interactions. 

Moreover, the pseudorapidity distributions of the 
anisotropy parameter ^(f?) of final charged hadrons has 
been found to be intimately related to the corresponding 



dynamics. Fast equilibration due to parton recombina- 
tion is necessary in order to obtain V2(i r ]) profiles which 
are clearly peaked at midrapidity as seen in the data. 
We want to note that the other observables studied with 
the standard QGSM [22j turned out to be essentially ro- 
bust against the inclusion of quark recombination. In 
particular the particle species dependence of V2(j>t) is 
still reproduced with the extended QGSM. The rapidity 
distribution for the final hadrons and their directed flow 
vi(rj) are even considerably improved. 

In conjunction with our aforementioned results for the 
local equilibration behavior, this is from the microscop- 
ical point of view a strong indication for the creation of 
a strongly interacting partonic medium in Au+Au colli- 
sions at RHIC that is thermally equilibrated on a very 
short time scale. 
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